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SUMMARY 


Inviscid  flow  in  the  region  of  the  leading  edge  of  curved  airfoils 
with  attached  shock  waves  is  investigated.  Tables  and  charts  are  pre¬ 
sented  for  determining  the  surface -pres sure  gradient  and  the  shock-wave 
curvature  in  supersonic  flow  of  an  ideal  diatomic  gas.  The  results  cover 
a  range  of  Mach  numbers  from  1.5  to  infinity  and  deflection  angles  from 
zero  up  to  those  approaching  shock  detachment.  Calculations  of  surface- 
pressure  gradient  and  shock-wave  curvature  are  also  made  for  curved  air¬ 
foils  in  supersonic  flow  of  a  calorically  imperfect,  diatomic  gas.  These 
calculations  are  quantitatively  applicable  in  cases  where  the  air  temper¬ 
atures,  downstream  of  the  shock  wave,  do  not  exceed  about  5 000°  Rankine . 

When  flow  conditions  approach  those  at  which  shock  waves  detach  from 
airfoils,  the  surface-pressure  gradient  and  shock-wave  curvature  vary 
widely  from  the  values  predicted  by  a  generalized  shock-expansion  method. 
Otherwise,  the  use  of  the  shock-expansion  method  introduces  only  small 
errors,  particularly  in  the  case  of  ideal  gas  flows.  The  effect  of 
caloric  imperfections  in  air  is  to  increase  these  errors. 

An  approximate  procedure  for  determining  the  flow  field  a  short 
distance  downstream  of  the  leading  edge  is  also  presented. 


INTRODUCTION 


Flow  in  the  region  of  the  leading  edge  of  an  airfoil  with  attached 
shock  wave  has  been  studied  by  numerous  investigators,  including  Crocco 
(reference  l),  and  more  recently  Munk  and  Prim  (reference  2),  Schaefer 
(reference  3)>  and  Thomas  (reference  4).  With  the  assumptions  that  the 
flow  is  steady,  two-dimensional,  and  inviscid,  and  that  air  behaves  as  an 
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ideal  (diatomic)  gas,  it  has  been  found  that  effects  on  the  flow  of 
interaction  between  the  shock  wave  emanating  from  the  leading  edge  and 
other  disturbances  emanating  from  the  surface  can  be  calculated.  It  is 
to  be  expected,  therefore,  that  within  the  limitations  of  these  assump¬ 
tions,  such  parameters  as  surface -pressure  gradient  and  shock-wave 
curvature  at  the  leading  edge  of  an  airfoil  can  be  accurately  predicted. 
The  first  purpose  of  the  present  report  is  to  present  values  of  these 
parameters  for  a  wide  range  of  Mach  numbers  and  flow  deflection  angles, 
thereby  supplementing  and  extending  the  results  presented  in  reference  2. 

Supersonic  flow  about  airfoils  was  investigated  more  generally  by 
Eggers  and  Syvertson  (reference  5)  an<^  it  was  found  that  caloric  imper¬ 
fections  in  air  influence  the  flow  appreciably  at  higher  Mach  numbers 
where  the  temperatures  of  the  disturbed  air  are  well  above  ambient 
temperatures.  They  used  a  generalized  method  of  characteristics  to  take 
into  account  these  imperfections  at  temperatures  up  to  about  5000° 
Rankine.  It  was  further  shown  that  the  shock-expansion  method,  general¬ 
ized  to  include  major  effects  of  interaction  between  shock  waves  and 
disturbances  incident  thereon  and  to  account  for  effects  of  caloric 
imperfections,  should  predict  the  flow  field  about  an  airfoil  with 
reasonable  accuracy.  The  advantage  of  this  method  over  the  method  of 
characteristics  is,  of  course,  its  relative  simplicity. 

Flow  in  the  region  of  the  leading  edge  of  an  airfoil  was  not  con¬ 
sidered  in  detail,  however,  in  reference  5-  An  additional  objective  of 
the  present  paper  is,  then,  to  determine  effects  of  caloric  imperfections 
on  flow  in  this  region,  using  a  method  analogous  to  that  of  Munk  and  Prim 
for  ideal  gas  flows,  particular  attention  being  given  to  surface-pressure 
gradient  and  shock-wave  curvature  at  the  leading  edge.  Predictions  of 
these  quantities  by  the  generalized  shock-expansion  method  are  also  pre¬ 
sented  for  both  ideal  and  calorically  imperfect  gas  flows. 

Since  it  is  useful  to  know  the  flow  conditions  a  short  distance 
downstream  of  the  leading  edge  (e.g.,  to  start  characteristics  solutions 
for  flow  about  airfoils),  an  approximate  method  of  determining  this  flow 
field  is  also  discussed.  This  method  depends  on  a  knowledge  of  the 
surface -pres sure  gradient  and  shock-wave  curvature  at  the  leading  edge. 


SYMBOLS 


Ci,  C2  characteristic  coordinates,  feet 
K  curvature,  feet  1 

M  Mach  number  (ratio  of  local  velocity  to  local  speed  of  sound) 
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ratio  of  static  pressure  downstream  of  the  shock  wave  to  free- 
stream  static  pressure 

absolute  temperature ,  degrees  Rankine 

distance  measured  from  leading  edge  along  airfoil  surface,  feet 

rectangular  coordinates,  feet 

Mach  angle  ^arc  sin  i  degrees 

angle  between  shock  wave  and  flow  direction  just  downstream  of 
the  shock  wave,  degrees 

ratio  of  specific  heat  at  constant  pressure  to  specific  heat  at 
constant  volume 

angle  between  flow  direction  just  downstream  of  the  shock  wave 
and  flow  direction  of  the  free  stream,  degrees 

molecular  vibrational  energy  constant,  degrees  Rankine 
(5500°  R  for  air) 

ratio  of  the  shock-wave  curvature  to  that  given  by  the  shock- 
expansion  method 

ratio  of  density  just  downstream  of  the  shock  wave  to  the  free- 
stream  density 

angle  between  flow  direction  of  free  stream  and  shock  wave, 
degrees 

ratio  of  surface -pres sure  gradient  to  that  given  by  the  shock- 
expansion  method 

Subscripts 

free -stream  conditions 

ideal  gas  quantities 

conditions  at  the  leading  edge  immediately  downstream  of  the 
shock  wave 


s  conditions  along  the  shock  wave 

w  conditions  along  the  airfoil  surface 
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ANALYSIS  OF  FLOW  IN  TEE  REGION  OF  THE  LEADING  EDGE 


Calculation  of  Surface -Pres sure  Gradient 
and  Shock-Wave  Curvature 


A  curved  two-dimensional  airfoil  surface  creates  pressure  disturb¬ 
ances  which  alter  the  inclination  of  an  attached  leading-edge  shock 
wave.  Figure  1  illustrates  the  geometry  of  this  phenomenon  for  a  convex 
airfoil.  The  pressure  disturbances  (expansion  waves  for  a  convex 
airfoil)  travel  along  first -family  Mach  lines  Ci  and  interact  with  the 
oblique  shock  wave.  In  addition  to  changing  the  inclination  of  the 
shock  wave,  this  interaction  produces  another  system  of  disturbances 
which  travel  along  second-family  Mach  lines  C2  from  the  shock  wave  to 
the  body. 


Method  of  characteristics.-  An  exact  solution  for  the  surface- 
pressure  gradient  and  shock-wave  curvature  at  the  leading  edge  may  be 
determined  by  the  method  of  characteristics  in  the  following  manner. 
It  is  clear,  referring  to  figure  1,  that  the  difference  in  flow  angle 
between  points  A  and  C  given  by  the  compatibility  equations  (see 
reference  2) 


and 


f  £6  \ 


sin  2  a 
2  7  P 

sin  2  a 
2  7  P 


(1) 

(2) 


along  the  path  ABC  must  equal  that  determined  by  the  airfoil  surface 
from  A  to  C.  Also,  the  difference  in  pressure  between  points  B  and  D 
given  by  these  compatibility  equations  along  the  path  BCD  must  equal 
that  determined  by  the  change  in  shock -wave  inclination  between  B  and  D. 
In  reference  2,  these  conditions  were  employed  at  the  leading  edge  to 
obtain  equations,  in  a  simple  parametric  form,  for  determining  the 
surface -pres sure  gradients  and  the  shock-wave  curvatures. 


These  equations  can  be  written  in  the  form 

1  dP  Z  +  1 


* 


* 


(3) 
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for  the  surface -pres sure  gradient,  and 

Js 

Kw 

for  the  shock-wave  curvature,  where 

dP/dCg  _  (d8/dP)CsX+  1  T  sin  (a  +  a  -  5)  1 

Z  ^P/BCi  (dP/d&)  (dS/dP)-"  -  1  sin  (a  -  a  +  5) 

L_  S  -2  —  — 


Z  sin  (a  -  0  +  &)  -  sin  (a  +  0  -  5) 


(4) 


(5) 


A  procedure  for  evaluating  equations  (3)  and  (4)  for  a  calorically 
imperfect,  diatomic  gas,  as  well  as  for  an  ideal  gas,  is  presented  in  the 
appendix  of  this  report.  The  values  of  surface -pres sure  gradient  and 
shock-wave  curvature  at  the  leading  edge  determined  using  these  equa¬ 
tions  are  exact  for  two-dimensional,  steady,  inviscid  flows. 

Two  applications  of  these  equations  will  he  considered:  First,  the 
accuracy  of  approximate  methods  of  calculating  the  flow  field  about  a 
curved  airfoil  can  be  evaluated  a.t  the  leading  edge  by  comparing  the 
values  of  surface -pressure  gradient  and  shock-wave  curvature  predicted 
by  these  approximate  methods  to  the  values  obtained  using  equations  (3) 
and  (4).  Second,  the  pressure  gradient  and  shock-wave  curvature  can  be 
used  to  determine  the  approximate  flow  field,  between  the  body  and  the 
shock  wave,  a  short  distance  downstream  of  the  leading  edge  and  thereby 
obtain  the  initial  points  of  a  characteristics  solution. 

Shock-expansion  method.-  An  approximate  solution  for  the  surface- 
pressure  gradient  and  shock-wave  curvature  may  be  obtained  using  the 
shock-expansion  method.  Eggers  and  Syvertson  (reference  5-)>  in  their 
general  discussion  of  the  entire  flow  field  about  airfoils,  found  that 
disturbances  incident  on  an  oblique  shock  wave  were  almost  entirely  con¬ 
sumed  in  changing  the  shock-wave  inclination,  and  that  the  reflected 
disturbances  were  weak  enough,  in  most  cases,  to  be  neglected.  Thus,  as 
a  first  approximation,  the  deflection  angle  and  pressure,  but  not  neces¬ 
sarily  the  Mach  number  and  entropy,  can  be  considered  constant  along 
first-family  Mach  lines.  The  use  of  this  assumption  leads  to  a  gener¬ 
alized  shock-expansion  method  for  calculating  the  entire  flow  field  about 
an  airfoil.  The  application  of  this  approximate  method  to  the  analysis 
of  flow  in  the  region  of  the  leading  edge  is  now  considered. 


If  the  disturbances  reflected  from  the  shock  wave  can  be  neglected. 


dP 
dC  1 
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Hence,  from  equation  (5),  Z  must  be  large  compared  to  unity.  If  this 
approximation  is  made,  the  expression  for  the  surface -pres sure  gradient 
(equation  (3))  reduces  to 


1  dP  1 

Kw  dW  rd&/dP)c2 

Similarly,  the  expression  for  the  shock-wave  curvature  (equation  (4)) 
reduces  to 


Ks  sin  (a  -  a  +  5) 
Kw 


dS 

da 


■5—  sin  a 


(8) 


It  should  be  realized  that  the  flow  field  is  determined  by  the 
basic  flow  equations  in  conjunction  with  the  shock  wave  and  airfoil 
surface  as  boundary  conditions.  Thus,  the  additional  requirement  for 
this  shock-expansion  method  of  zero  pressure  gradient  along  first-family 
Mach  lines  means  that  one  of  the  flow  relations  cannot  be  satisfied 
exactly  (i.e.,  the  flow  field  is  overdetermined).  Equations  (7)  and  (8) 
satisfy  the  shock  relations  and  the  airfoil  surface  as  boundary  condi¬ 
tions;  however,  the  compatibility  equations  are  only  approximately 
satisfied. 


The  error  in  surface -pressure  gradient  associated  with  neglecting 
the  reflected  disturbances  might  be  expected  to  be  largest  in  the  region 
of  the  leading  edge  of  a  curved  airfoil  due  to  the  close  proximity  of 
the  shock  wave  and  the  surface.  The  magnitude  of  the  error  in  this  region 
may  be  deduced,  of  course,  from  the  ratios  of  values  of  surface -pres sure 
gradient  and  shock-wave  curvature  given  by  the  characteristics  method  to 
those  given  by  the  shock-expansion  method.  The  surface-pressure-gradient 
ratio  and  shock-wave -curvature  ratio  can  be  written  (using  equations  (3) 
and  (7)) 


\) 1  _  z  +  1 
Z  -  1 

and  (using  equations  (A)  and  (8)) 


Z  sin  (a  -  q  +  5)  -  sin  (a+g-S) 
(z  -  l)  sin  (a  -  a  +  5) 


(9) 


(10) 


respectively.  A  procedure  for  evaluating  equations  (9)  and  (10)  for 
flow  of  a  calorically  imperfect  gas,  as  well  as  for  an  ideal  gas,  is 
presented  in  the  appendix  of  this  report.  (The  application  of  equa¬ 
tion  (9)  for  ideal  gas  flow  has  already  been  given  In  reference  2.) 


MCA  TN  2729 


7 


Calculation'  of  the  Flow  Field  a  Short  Distance 
Downstream  of  the  Leading  Edge 


A  knowledge  of  the  values  of  flow  parameters  along  a  line  from  the 
"body  to  the  shock  wave  is  useful  in  starting  characteristics  or  shock- 
expansion  solutions  for  the  entire  flow  field.  A  method  is  now  pre¬ 
sented  for  approximating  these  parameters  along  a  first-family  Mach  line 
a  short  distance  downstream  of  the  leading  edge.  In  the  region  of  the 
leading  edge,  the  airfoil,  the  shock  wave,  and  a  first-family  Mach  line 
are  approximated  by  circular  arcs,  as  shown  in  figure  2.  The  equations 
of  the  airfoil  surface  and  the  shock  wave  can  then  be  written  (knowing 
the  slope  and  curvature  at  the  leading  edge) 


cos 


Yw  =  “ 


SN  J1  -  (kwn  xw  ~  sin  5^ 


K 


-wN 


K, 


(ID 


and 


Ys  =  - 


cos 


ctn 


K, 


SN 


/  -  (ksn  *s  -  sin  oN)‘ 
_ 
sN 


(12) 


respectively.  If  it  is  assumed  that  the  pressure  gradient  is  constant 
along  the  airfoil  surface  for  a  short  distance  downstream  of  the  leading 
edge,  the  pressure  distribution  and.  thus  all  other  flow  conditions  at 
the  surface  are  known  since  the  flow  conditions  at  the  leading  edge  are 
known  (see  previous  analysis).  Since  the  shock  wave  is  approximated  by 
a  circular  arc,  the  flow  conditions  just  downstream  of  the  shock  wave  are 
also  known.  Thus,  the  flow  parameters,  and  consequently  the  slopes  of  a 
first-family  Mach  line,  will  be  known  at  its  intersections  with  the  body 
and  the  shock  wave.  The  problem,  then,  is  to  determine  the  curvature, 
and  thus  the  equation  of  this  Mach  line.  The  following  analysis  is 
devoted  to  this  matter. 


To  the  accuracy  of  the  present  analysis,  the  flow  angle  5  and  the 
Mach  angle  a  vary  linearly  with  distance  along  the  Mach  line.  Thus 
(see  fig.  2) 


yB  ~  yA 
XB  “  XA 


(aw  +  ^w)a  +  (as  +  $s) 


B 


(13) 
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where  the  local  flow  angle  and  local  Mach  angle  at  the  airfoil  surface 
can  he  written  (as  a  first  approximation) 


c  ~  c  ,  (  d&  )  ( dw'') 

5a  -  5W  +  Vdw/  i 

N  v  XN 


XA  =  SN  +  kwn  xA  sec  % 


„  ^  f  da 

aA  =  aN  +  (  - 


'  i  dP\ 

dW  )  KwN  xA  sec  5N 
'  N 


respectively. 


Similarly,  at  the  shock  wave, 

6b  ~  aH  +  (f) 


Ksn  xB  sec  °W 


cy  ,  /  da 

^  W  Vdc 


KSn  xb  Bec  aw 


I'i 

It  is  possible  to  determine  ,  (ff )  ,  and  for 

SN  SN  VW 

imperfect  gas  flow  as  well  as  for  ideal  gas  flow.  Thus,  an  implicit 
relation  between  xs  and  xw  can  be  found  from  equations  (ll)  through 

The  equation  of  the  circular  arc  approximating  the  first-family 
Mach  line  may  be  written,  for  any  initial  xw  (knowing  the  slopes  at  the 
intersection  of  this  line  with  the  shock  wave  and  the  body) 

„  _  c°s  (Kww  ^  “  sin 

"  v.  +  P - -  - 


« 
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At  point  B,  the  intersection  of  the  shock  wave  and  the  first-family  Mach 
line,  equation  (l8)  may  he  solved  for  Kg  ,  using  the  implicit  relation 

between  xs  and  xw.  Substituting  this  value  of  Kq.,  into  equation  (18), 
an  equation  relating  yQ^  and  x^.  along  the  first -family  Mach  line  can 

be  found  for  every  value  of  xw. 

The  pressure  and  flow  angle  at  A  and  B  and  the  changes  in  pressure 
and  flow  angle  from  A  to  B  may  be  determined.  The  change  in  pressure 
with  change  in. flow  angle  should  agree  (to  the  desired  degree  of  accu- 
racy)  with  that  calculated  using  the  basic  compatibility  relation, 
namely. 


B 

2-  _  ?(PA  +  PB) 
sin  (aA  +  aB) 

A 

Aay  significant  difference  should  be  diminished  by  choosing  point  A 
closer  to  the  leading  edge. 

«  The  angles  8  and  a  are  known  at  each  point  along  this  first- 

family  Mach  line,  and  the  pressure  variation  may  be  determined  using 
the  characteristics  compatibility  equation  (equation  (l)).  Hence  the 
entropy  variation  can  also  be  determined.  Therefore,  all  necessary 
flow  parameters  are  defined  along  this  line,  rendering  it  a  good  start¬ 
ing  line  for  a  characteristics  solution  for  the  entire  flow  field. 


AP 

AS 


DISCUSSION 
Ideal  Gas  Flows 


Surface-pressure  gradient.-  The  results  of  the  calculations  (using 
equation  (A4)  of  the  appendix)  of  surface -pressure  gradient  are'  presented 
in  table  I  and  figures  3  and  4.  The  values  presented  in  table  I  are  for 
a  range  of  Mach  numbers  from  1.5  to  °°  and  for  leading -edge  deflection 
angles  from  0°  to  4-5°.  Where  no  value  appears  in  the  table,  the  flow 
behind  the  shock  wave  is  subsonic.  Additional  calculations,  not  pre¬ 
sented  in  the  table,  were  made  to  permit  the  accurate  plotting  of  curves 
to  the  point  where  the  flow  behind  the  shock  wave  attains  sonic  velocity.1 

Charts  were  also  presented  for  surf ace -pressure  gradient,  surface- 
pressure  -gradient  ratio,  and  shock-wave  curvature  for  ideal  gas  flows 
in  reference  2;  however,  the  results  given  in  the  present  report  are 
_  somewhat  more  extensive . 
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Surface-pressure-gradient  ratio.-  The  results  of  the  calculations 
of  surface-pressure-gradient  ratio  (equation  (A6))  are  presented  in 
table  II  and  figures  5  and  6  in  a  manner  similar  to  the  presentation  of 
the  surface -pres sure  gradient.  From  these  results,  it  is  seen  that, 
except  near  detachment,  the  surf ace -pres sure -gradient  ratio  varies 
between  approximately  0.88  and  1.12  for  all  Mach  numbers  including  <». 
Below  a  Mach  number  of  4  the  ratio  varies  only  from  O.98  to  1.02, 
except  near  detachment.  Therefore,  little  error  will  result  from  the 
use  of  the  shock-expansion  method  for  the  surface -pres sure  gradient  at 
the  leading  edge  of  airfoils  for  Mach  numbers  less  than  4.  Near  detach¬ 
ment,  for  Mach  numbers  greater  than  approximately  2,  the  surface- 
pressure  -gradient  ratio  attains  a  very  large  range  of  values.  The  maxi¬ 
mum  value  of  the  ratio  increases  with  Mach  number.  The  use  of  the 
surface -pres sure  gradient  at  the  leading  edge  as  given  by  the  shock- 
expansion  method  would  result  in  appreciable  error  near  detachment 
conditions. 

It  should  be  noted  that,  for  the  particular  case  of  infinite  free- 
stream  Mach  number  and  zero  deflection  angle,  the  pres sure -gradient 
ratio  ^  is  double -valued.  From  equation  (A6),  it  is  apparent  that  + 
is  1  for  zero  deflection.  Yet,  at  infinite  free-stream  Mach  number, 

\|r  approaches  0.8822  as  the  deflection  angle  approaches  zero. 

Tile  flow  along  the  surface  is  isentropic.  Hence,  it  can  be  shown 
that  i|r,  the  ratio  of  surface -pressure  gradient  to  that  given  by  the 
generalized  shock-expansion  method,  is  also  the  velocity-gradient  ratio, 
the  Mach  number  gradient  ratio,  the  Mach  angle  gradient  ratio,  the 
density-gradient  ratio,  and  the  temperature -gradient  ratio.  Any  of  these 
gradients  may  be  found,  then,  by  calculating  the  gradient  using  the 
shock-expansion  method  and  applying  the  appropriate  value  of  ty.  This 
property  of  the  ratio  makes  it  useful  in  the  application  of  the 
method  of  characteristics  with  any  of  the  coordinate  systems  commonly 
employed  in  the  compatibility  equations. 

Shock-wave  curvature.-  The  results  of  the  calculations  of  shock- 
wave  curvature  (equation  (A5))  are  presented  in  table  III  and  figures  7 
and  8.  It  is  clear  from  equation  (A5)  that  Kb/kw  is  zero  for  zero 
deflection  angle.  However,  when  the  free-stream  Mach  number  is  infinite, 
Ks/Kw  approaches  0.8003  as  the  deflection  angle  approaches  zero.  Thus, 
Ks/Kw  is  double -valued  at  M  =  00,  6  =  0. 


eThe  precise  value  of 
3The  limiting  value  of 


t  for  this  limit  is  2  y  -1 

K  /K.  is  +  1)-g 
Ks/Kw  is  4(27  -1)* 
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Shock-wave -curvature  ratio . -  The  shock-wave -curvature  ratio  for 
various  free-stream  Mach  numbers  and  deflection  angles  has  been  calcu¬ 
lated  by  means  of  equations  (8)  and  (A5) .  The  results  of  these  calcula¬ 
tions  are  presented  in  table  IV  and  figures  9  and  10.  It  is  interesting 
to  note  that  the  curves  in  figure  9  are  similar  to  the  reciprocals  of  the 
curves  shown  in  figure  5*  Indeed,  this  conclusion  can  be  deduced  by 
comparing  equations  (9)  and  (10).  Except  near  detachment,  the  curvature 
ratio  varies  from  0.92  to  1.08  for  all  Mach  numbers  including 
Thus,  only  moderate  errors  would  result  from  using  the  value  of  shock- 
wave  curvature  given  by  the  shock-expansion  method  for  all  flow  condi¬ 
tions  except  near  detachment.  For  zero  deflection  angle,  it  is  apparent 
from  equation  (10)  that  k  must  be  1.  However,  for  infinite  free-stream 
Mach  number,  k  approaches  1.0724  as  the  deflection  angle  approaches 
zero.  Thus  k  is  also  double -valued. 


Calorically  Imperfect  Gas  Flows 

The  temperature  ratio  across  an  oblique  shock  wave  increases  with 
increasing  Mach  number  and  leading-edge  slope,  as  shown  in  figure  11. 

As  the  temperature  behind  the  shock  wave  increases,  the  behavior  of  the 
air  diverges  from  that  of  an  ideal  gas.  Below  800°  Rankine,  the  diver¬ 
gence  is  not  significant,  and  the  equations  for  ideal  gas  flow  can  be 
applied  with  only  minute  errors  resulting. 

Above  8 00°  Rankine,  the  energy  of  the  vibrational  degrees  of  freedom 
of  the  gas  molecules  is  appreciable  and  becomes  greater  with  increasing 
temperature .  As  a  result,  the  ratio  of  specific  heats  varies  with  tem¬ 
perature,  as  discussed  by  Eggers  in  reference  6. 

As  mentioned  previously,  equations  were  developed  in'  reference  5 
permitting  the  calculation  of  the  flow  of  a  diatomic,  calorically 
imperfect  gas  about  two-dimensional  airfoils.  These  equations  apply  for 
temperatures  up  to  the  order  of  5000°  Rankine  and  form  the  basis  for  the 
extension  of  the  solution  for  surface -pres sure  gradient,  surface -pressure 
gradient  ratio,  shock-wave  curvature,  and  shock-wave -curvature  ratio  to 
the  case  of  calorically  imperfect  gases.  For  a  free-stream  temperature 
,  of  500°  Rankine,  the  shaded  area  between  lines  of  constant  temperature 
ratio  of  1.5  and  10,  in  figure  11,  represents  the  approximate  range  of 
conditions  for  which  the  method  developed  in  this  report  for  the  flow  of 
a  diatomic,  calorically  imperfect  gas  would  apply.  The  range  shown  in 
figure  11  is  only  approximate,  since  it  was  calculated  using  ideal  gas- 
flow  equations. 


4 The  exact  limiting  value  of  k  is 


7  +  1 


2(27  -  1)  (l  -  S' 


\ 
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The  excitation  of  the  vibrational  degrees  of  freedom  of  the  gas 
molecules  requires  a  finite  number  of  collisions,  causing  the  well-known 
heat-capacity  lag  discussed  by  Bethe  and  Teller  in  reference  7-  The 
flow  distance  (i.e.,  along  the  streamline)  required  to  establish  equi¬ 
librium  conditions  is  usually  small  in  dense  air  and  will  be  considered 
infinitesimal  in  this  report.5 

Since  the  free -stream  static  temperature  is  an  additional  parameter 
in  calculations  of  flow  of  imperfect  gases,  only  a  limited  number  of 
calculations  of  (l/Kv)  (dP/dW) ,  \|r,  Ks/kw,  and  k  were  made.  The  purpose 
of  these  calculations  is  to  compare  the  variations  of  these  quantities 
with  the  values  as  given  by  the  ideal -gas -flow  computations.  The  calcu¬ 
lations  followed  the  procedure  described  in  the  appendix.  A  free -stream 
static  temperature  of  500°  Rankine  was  used.  The  results  of  these  calcu¬ 
lations  are  presented  in  table  V  for  various  Mach  numbers  and  leading - 
edge  deflection  angles. 

The  surface-pressure  gradients  for  an  ideal  gas  and  for  a  calorically 
imperfect,  diatomic  gas  are  compared  in  figure  12.  In  all  cases  calcu¬ 
lated,  the  gradient  is  smaller  for  the  imperfect  gas  and  diverges  gradu¬ 
ally,  with  increasing  free-stream  Mach  number  and  deflection  angle,  from 
the  value  of  the  gradient  for  an  ideal  gas.  This  divergence  is  consist¬ 
ent  with  the  increasing  effects  of  the  caloric  imperfections  due  to  the 
increasing  temperature  behind  the  shock  wave. 

The  surf ace -pres sure -gradient  ratio  for  the  imperfect  gas  is 
compared  in  figure  13  with  the  ratio  for  an  ideal  gas.  The  angle  (3 
between  the  shock  wave  and  the  leading  edge  of  the  airfoil  is  smaller 
for  imperfect  gas  flows  than  for  ideal  gas  flows,  and  shock-wave 
expans ion -wave  interaction  effects  are  increased.  As  a  result,  the 
surface-pressure-gradient  ratio  is  smaller  for  imperfect  gas  flows  than 
for  ideal  gas  flows.  This  difference  in  the  angle  P  is  a  function  of 
the  temperature  behind  the  shock  wave.  Thus,  the  surface-pressure- 
gradient  ratio  for  imperfect  gas  flows  diverges  from  the  ratio  for  ideal 
gas  flows  with  increasing  Mach  number  and  deflection  angle. 

A  divergence  with  Mach  number  and  deflection  angle  is  also  apparent 
in  figure  14  in  which  is  compared  the  shock-wave  curvature  for  an  ideal 
gas  and  a  calorically  imperfect,  diatomic  gas.  This  divergence  is 
compatible  with  the  change  in  surface -pres sure  gradient  due  to  the  caloric 
imperfections  of  the  gas. 


5 The  dissociation  of  air  at  temperatures  greater  than  about  5000°  Rankine 
is  also  discussed  in  reference  7-  This  phenomenon  is  not  considered  in 
the  present  report. 
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The  shock-wave -curvature  ratio  for  a  calorically  imperfect, 
diatomic  gas  and  this  ratio  for  an  ideal  gas  are  shown  in  figure  15. 
Again  it  is  seen  that  the  effect  of  caloric  imperfection  is  to  increase 
the  effects  of  shock-wave  expansion-wave  interaction. 


CONCLUDING  REMARKS 


The  supersonic  flow  field  in  the  region  of  the  leading  edge  of 
curved  two-dimensional  airfoils  has  been  investigated. 

A  set  of  tables  and  charts  has  been  presented  giving  the  surface- 
pressure  gradient  and  shock-wave  curvature  at  the  leading  edge  for  an 
ideal  gas  for  a  large  range  of  free -stream  Mach  numbers  and  leading-edge 
deflection  angles.  Previously  derived  equations  for  surface -pres sure 
gradient  and  shock-wave  curvature  have  been  extended  to  permit  the 
determination  of  these  quantities  for  flow  of  a  calorically  imperfect, 
diatomic  gas.  Both  the  previously  derived  equations  and  the  extended 
equations  have  been  modified  so  that  the  shock-wave  curvature  and 
surface -pres sure  gradient  are  consistent  with  the  assumptions  of  a 
generalized  shock-expansion  method.  Values  of  surface -pres sure  gradient 
and  shock-wave  curvature  have  been  compared  with  the  values  given  by  the 
modified  equations.  It  has  been  shown  that,  except  near  conditions  at 
which  the  shock  wave  detaches  from  the  airfoil,  only  moderate  errors 
would  result  from  using  the  values  given  by  shock-expansion  equations 
for  flow  of  an  ideal  gas.  It  has  also  been  shown  that  caloric  imper¬ 
fections  of  air  increase  the  error. 

An  approximate  procedure  has  been  presented  for  determining  the 
flow  field  a  short  distance  downstream  of  the  leading  edge.  This  pro¬ 
cedure  depends  on  a  knowledge  of  the  surface -pres sure  gradient  and  shock- 
wave  curvature  at  the  leading  edge. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Mar.  31,  1952. 
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APPENDIX 


EVALUATION  OF  ( l/Kw) ( dP/dW) ,  Ks/Kw,  *,  and  K 


For  an  ideal  gas,  the  Mach  number,  Mach  angle,  shock-wave  angle, 

and  pressure  ratio  can  be  calculated  at  the  leading  edge  using  the 

standard  Rankine-Hugoniot  shock  relations  and  utilizing  the  free-stream 

Mach  number  and  leading-edge  deflection  angle.  Hence  (d5/dP)p  can  be 

l2 

determined  from  equation  (2)  of  the  analysis.  With  these  flow  parameters 
known,  the  only  terms  still  to  be  determined  in  equations  (3),  (^),  (5), 
(9),  and  (10)  are  (dP/d8)g  and  (dS/da)s. 


Now  (see  reference  2) 


f dp\  =  (dP/dg)s 
VdS/s  (d&/da ) s 


(Al) 


where 


(a.  • 


2  7  Mq  sin  2  o 

7  +  1 


(A2) 


and 


d5 

da 


2^y  (Mo2  sin  2  o)S  +  Mo2  sin2  ^  cos  2 


Mo2  sin  2  * 


f 


^  =  1  - 


1  + 


r  2  2  .  2  _ ">2 

+  Mq  sin  a 


-y7-1.-  M02  sin  2  a 

2(7-1)  0 


(A3) 


Equations  (A2)  and  (A3)  are  the  derivatives  of  the  forms  of  the  Rankine- 
Hugoniot  shock  relations  given  in  reference  2. 
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Using  the  standard,  forms  of  the  Rankine— Hugoniot  shock  relations, 
it  is  possible  to  transform  equations  (3),  (4),  and  (9)  (given  in  the 
analysis)  into 


_L  dP 
Kv  dw 


_  7  tan  B 
sin2  a 


1  +  l 


cos2  a 
cos2  p 


Mr 


sinc 


l) 


cos2  a 


tan2  p  1 

tan2  a  2  Vcos2  p 


M02  sin2  a 


2  7  Mq2  sin2  a  -(7-1) 


7  +  1 


(Ah) 


for  the  surface-pressure  gradient. 


^s 

Kw 


7  +  1 
V  cos  p 


1  - 


tan2  p 
tan2 


a 


tan2  p 
tan2  a 


+  2 


:0s2 


a 


cos2  p 


Mq2  sin2  a 


(A5) 


for  the  shock-wave  curvature,  and 


f  = 


1  /  cos  a 

1  +  —  I  - - —  + 


2  \cos2  P  Mq2  Sin2 


tan2  p  1 

- 5 —  +  — 

tan*  a  2 


cos2  a 
cos2  p 


+ 


Mq2  sin2  a 


tan  p 
tan  a 


(A 6) 


for  the  surface-pressure -gradient  ratio.  These  equations  are  similar  to 
those  given  by  Schaefer  in  reference  3  and  require  less  work  to  compute 
than  equations  (3),  (4),  and  (9). 


For  a  calorically  imperfect,  diatomic  gas,  the  standard  shock 
relations  obviously  are  not  applicable.  However,  expressions  for  the 
relations  for  flow  of  a  calorically  imperfect  gas  through  an  oblique 
shock  wave  have  been  developed  in  reference  5* 
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Thus  we  have,  as  a  function  of  the  local  static  temperature  and 
free-stream  conditions  _  _ 


where  7^  =  1.400  and  6  =  5500°  R* 


The  Mach  number  immediately  downstream  of  the  shock  can  be  written 


Zo^of  1  n  (± 

2  7i“l  \ 


±( — A _ L_\ 

To(  9/T0  9/T3  J 

\  e  -1  e  -V 


(A8) 


The  relations  for  the  pressure  and  density  ratios  across  the  shock  are 


Ps  = 


(1+7sM82)-  ^ 
J-s 


(l+70Mo  )  +  V 


r 


(1+7sMb2)-  ^  (1+7QMQ2) 


4_Tc 

To 


(A9) 


and 


_1_  Ts 
Ps  T0ps 


(A10) 


respectively.  The  expressions 
deflection  angle  are 


sin2  0 


and 


for  the  shock-wave  angle  and  the  flow 


^  J® 

y o  To  \m0/ 


(All) 


cot  6S 


tan  a 


(A12) 
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respectively.  Equations  (A7)  through  (A12)  are  valid  for  temperatures 
downstream  of  the  shock  wave  up  to  the  order  of  5000°  R. 

The  primary  variable  is  the  temperature,  immediately  downstream  of 
the  shock  wave.  Hence,  we  can  write 


(cLP/dS)s 


(ap/aT)s 

(dS/dT)s 


(A13) 


and 


(dS/d0)s 


(dS/dT)s 

(dcj/dT)s 


(All*) 


Differentiation  of  equation  (A9)  yields 
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8. 


9- 


10. 


11. 


12. 


13. 


14. 


15- 


16. 


17. 


18. 


19. 


20. 


21. 


Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 

Calculate 


(equation  (A17)) 
(equation  (Al6) ) 
(equation  (A15) ) 

(equation  (A20)) 
(equation  (A19)) 
(equation  (Al8)) 
(equation  (A13)) 
(equation  (2)) 
(equation  (5)) 
(equation  (3)) 
(equation  (9)) 
(equation  (Al4)) 
(equation  (4)) 
(equation  (10)) 


t 


MCA  TN  2729 


21 


REFERENCES 


1.  Crocco,  Luigi;  Singolarita  della  Corrente  Gassosa  Iperacustica 

Nell’Intorno  di  una  Prora  a  Diedro.  L'Aerotecnica,  vol.  17,  no.  6, 
June  1937. 

2.  Munk,  M.  M. ,  and  Prim,  R.  C.:  Surface -Pres sure  Gradient  and  Shock 

Front  Curvature  at  the  Edge  of  a  Plane  Ogive  with  Attached  Shock 
Front.  Jour.  Inst.  Aero.  Sci.,  vol.  15,  no.  11,  Nov.  1948, 
pp.  691-695. 

3.  Schaefer,  M.:  The  Relation  Between  Wall  Curvature  and  Shock  Front 

Curvature  in  Two  Dimensional  Gas  Flow.  A.F.,  AMC,  Tech.  Rept.- 
F-TS-1202-IA, (GDAM  A9-T-9)  Brown  University.  (Technische 
Hochschule  Dresden,  1942.  Peenemunde  Archiv  44/8.) 

4.  Thomas,  T.  Y.:  Calculation  of  the  Curvature  of  Attached  Shock 

Waves.  Jour.  Math,  and  Physics,  vol.  27,  no.  4,  Jan.  1949, 

pp.  .279-297. 

5*  Eggers,  A.  J.,  Jr.,  and  Syvertson,  Clarence  A.:  Inviscid  Flow  About 
Airfoils  at  High  Supersonic  Speeds.  NACA  TN  2646,  1952. 

6.  Eggers,  A.  J.,  Jr.:  One -Dimensional  Flows  of  an  Imperfect  Diatomic 

Gas.  NACA  Rep.  959,  1950.  (Formerly  NACA  TN  l86l,  1949) 

7.  Bethe,  H.  A.,  and  Teller,  E.:  Deviations  from  Thermal  Equilibrium 

in  Shock  Waves.  Aberdeen  Proving  Ground,  Maryland.  Ballistic 
Research  Lab.  Rep.  No.  X-117,  1945. 


TABLE  II.-  SURFACE -PRESSURE -GRADIENT  RATIO 


NACA 


TABLE  III.-  SHOCK-WAVE  CURVATURE 


2k 


NACA  TN  2729 


TABLE  IV.-  SHOCK-WAVE -CURVATURE  RATIO 
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TABLE  V.-  SURFACE -PRESSURE  GRADIENT,  SURFACE-PRESSURE-GRADIENT 
RATIO,  SHOCK-WAVE  CURVATURE,  AND  SHOCK-WAVE -CURVATURE 
RATIO  FOR  A  CALORICALLY  IMPERFECT,  DIATOMIC  GAS 
[T0  =  500°  RANKINE] 


Mo 

S° 

VT0 

1  dP 
Kw  dW 

* 

KsAw 

K 

3 

30.4757 

2.012 

18.36 

1.0262 

1.0212 

0.9828 

5 

20.4748 

2.138 

30.94 

.9675 

.5866 

1.0312 

5 

31.0968 

3.127 

42.63 

.9548 

.7987 

1.0313 

5 

41.0064 

4.300 

60.76 

1.0419 

1.5991 

1.0000 

10 

10.0129 

2.108 

63.87 

.9622 

.5266 

1.0394 

10 

20.5239 

4.400 

123.6 

.9111 

.7137 

1.0621 

10 

31.2637 

7.640 

169.8 

.8889 

.8363 

1.0636 

10 

42.9772 

a12.000 

205.2 

.9178 

1.1678 

1.0365 

20 

10.6610 

4.700 

274.1 

.9024 

.6912 

1.0685 

20 

20.9696 

&12 . 600 

515.0 

.8592 

.8o4i 

1.0787 

20 

32.5087 

a25.8oo 

705-1 

.8492 

.8915 

1.0753 

20 

39.1836 

a35.000 

770.4 

.8559 

.9858 

1.0664 

aVaTues  of  temperature  greater  than  5000°  R  downstream 
of  the  shock  wave.  '\nacaD-'7’ 


wave  approximated  by  circular  arcs. 
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Leading-edge  deflection  angle,  8,  degrees 
(a)  Mo,  1.5  to  6 


Figure  3.— Variation  of  surface-pressure  gradient  with  leading-edge 
deflection  angle  for  various  free-stream  Mach  numbers. 
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Figure  3.— Concluded. 
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Figure  4 . — Variation  of  surface-pressure  gradient  with 
free-stream  Mach  number  for  various  leading-edge 
deflection  angles. 
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Figure  4. — Concluded. 
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Figure  5-  Variation  of  surface-pressure-gradient  ratio  with 
deflection  angle  for  various  free-stream  Mach  numbers. 
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Free-stream  Mach  number.  Mo 
(a)  8,  0°  to  25° 


Figure  6r  Variation  of  surface-pressure-gradient  ratio  with  free-stream 
Mach  number  for  various  leading- edge  deflection  angles. 
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Figure  6r  Concluded. 


Shock-wave  curvature. 
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Leading-edge  deflection  angle,  8,  degrees 


Figure  7.  —  Variation  of  shock-wave  curvature  with  leading-edge 
deflection  angle  for  various  free- stream  Mach  numbers. 
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10-Variation  of  shock-wave-curvature  ratio  with  free-stream  Mach  number  for  various 

leading-edge  deflection  angles. 
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Free-stream  Mach  number,  Mo 


Figure  II. — Variation  of  leading-edge  deflection  angle  with  free-stream 
Mach  number  for  various  values  of  the  temperature  ratio. 
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Comparison  of  the  variation  of  surface-pressure  gradient  with  leading- edge 
a  angle  for  various  free- stream  Mach  numbers  for  an  ideal  gas  and  for  a  calorically 
t,  diatomic  gas. 
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Figure  14-  Comparison  of  the  variation  of  the  shock-wave  curvature  with  the  leading-edge 

deflection  angle  for  various  free-stream  Mach  numbers  for  an  ideal  gas  and  for  a 
calorically  imperfect,  diatomic  gas. 
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Figure  15.  - Comparison  of  the  variation  of  shock-wave-curvature  ratio  with  leading-edge 
deflection  angle  for  various  free- stream  Mach  numbers  for  an  ideal  gas  and  for  a 
calorically  imperfect \  diatomic  gas. 
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